Deformations and dilations of chaotic billiards, dissipation rate, 
and quasi-orthogonality of the boundary wavefunctions 
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We consider chaotic billiards in d dimensions, and study the matrix elements M„ m corresponding 
to general deformations of the boundary. We analyze the dependence of |M„ m | 2 on u) = (E n —E m )/h 
using semiclassical considerations. This relates to an estimate of the energy dissipation rate when the 
deformation is periodic at frequency ui. We show that for dilations and translations of the boundary, 
|Mn m | 2 vanishes like lj 4 as ui— >0, for rotations like ui 2 , whereas for generic deformations it goes to 
a constant. Such special cases lead to quasi-orthogonality of the eigenstates on the boundary. 



Chaotic cavities (billiards) in d dimensions are proto- 
type systems for the study of classical chaos and its fin- 
gerprints on the properties of the quantum-mechanical 
eigenstates. As the properties of static billiards are be- 
ginning to be understood, questions naturally arise about 
deformations and their time dependence. It is perhaps 
not widely appreciated that certain deformations are very 
special, and that there is a close connection between the 
quantum and classical mechanics of such deformations in 
the case of ergodic systems. In this paper, which takes 
a fresh approach to these issues, we explore a special 
class of deformations which do not 'heat' in the limit of 
small frequencies. We also establish a rather surprising 
relationship to a very successful numerical technique for 
finding billiard eigenfunctions. 

We start with the one-particle Hamiltonian Ho(r, p) = 
p 2 /2m + V(r), where m is the particle mass, r is the po- 
sition of the particle inside the cavity and p is the conju- 
gate momentum. We will take the limit V(r) —> oo out- 
side the cavity, zero otherwise, corresponding to Dirich- 
let boundary conditions. In this limit, the Hamiltonian 
is completely defined by the boundary shape. The vol- 
ume of the cavity we call V = L d . Upon quantization 
a second length scale A B = 2ir/k appears, where fc is 
the wavenumber. For simple geometries the typical time 
between collisions with the walls is r col ~ L/i>, where v 
is the particle speed. The energy is E = ^mv 2 . Upon 
quantization the eigenenergies are E n — (hk n ) 2 /2m. 
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FIG. 1. The two-dimensional (2D) billiards under numerical 
study. Left: Bunimovich quarter-stadium (a/R = 1). Right: Gen- 
eralized Sinai billiard (a/b = 2, Si = 0.2, 9 2 = 0.5). 

A powerful tool for the classical analysis is known as 
the 'Poincare section'. Rather than following trajectories 
in the full (r, p) phase-space, it is much more efficient 



to record only successive collisions with the boundary. 
This way we can deal with a canonical transformation 
(map) which is defined on a 2(d— 1) dimensional phase 
space. A similar idea is used in quantum-mechanics: By 
Green's theorem is is clear that all the information about 
an eigenstate ip(r) is contained in the boundary normal 
derivative function ip(s) = n-Vf/ 1 , where s is a (d—1) 
dimensional coordinate on the boundary, and n(s) the 
outward unit normal vector. 

However, unlike the classical case, the reduction to 
the boundary is not satisfactory. One cannot define 
an associated Hilbert space that consists of the bound- 
ary functions. In particular, the orthogonality relation 
{ipn\ipm) — &nm does not have an exact analog on the 
boundary. Still, the boundary functions 'live' in an effec- 
tive Hilbert space of dimension ~ (L/A B ) and it has 
been realized jl| that the following quasi-orthogonality 
relation holds. Define an inner product 

M nm = — ^ I <p n (s)(p m (s) (n-D)ds (1) 

where D(s) = r(s) is the displacement field correspond- 
ing to dilation (about an arbitrary origin), and k n ps 
k m ps k [p| . It is well known that the normalization con- 
dition (4>n\ipn) = 1 implies M nn = 1. We give a proof of 
this exact result in the Appendix. On the other hand the 
off-diagonal elements are only approximately zero || . 

The main purpose of this Letter is to study the band 
profile of the matrix M nm for a general displacement field 
D(s). In particular we want to understand why for spe- 
cial choices of D(s), notably dilations, we have quasi- 
orthogonality. Later we will explain that M nm can be 
interpreted as the matrix element of a perturbation STi. 
associated with a deformation of the boundary, such that 
(n-D)Jx is the normal displacement of a wall element, 
given a control parameter Sx. In the following two para- 
graphs we explain the main motivations for our study. 

The matrix elements \M nm \ 2 determine the rate of 
irreversible energy absorption by the particle {i.e. dis- 
sipation) due to external driving. Here 'external driv- 
ing' means time-dependent deformation of the boundary. 
Having exceptionally small |M„ m | 2 for special choices of 
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D(s), such as dilations, translations and rotations, im- 
plies exceptionally small dissipation rate ('non-heating' 
effect). This observation goes against the naive kinetic 
picture that the rate of heating should not depend on how 
we 'shake' the boundary. The special nature of transla- 
tions and rotations for uj = has been recognized in the 
context of nuclear dissipation . Our present approach 
allows us to analyze the non-heating effect present for di- 
lations as well, and provide the form of the low-frequency 
response of the system in all three cases (dilations, trans- 
lations and rotations). 

There is another good motivation to study this issue. 
Recently, a powerful technique for finding clusters of bil- 
liard eigenstates and eigenenergies has been found by 
Vergini and Saraceno with a speed typically ~ 10 3 
greater than previous methods. This efficiency relies on 
the above quasi-orthogonality relation, the associated nu- 
merical error being given by the deviation of M nm from 
S nm ■ Those authors tried to establish quasi-orthogonality 
using the identity M nm = S nm + [(fc^ - k 2 ,)/2k 2 ]B nm , 
with B nrn = (ijj n \r - V \ip m }, and by assuming Q that 
|-Bnm| ~ 0(1). However, a naive random wave argument 
would predict \B nm \ ~ 0(1/ X^' 1 ^ 2 . 

Fig. 2 displays the band profile \M nm \ 2 for three 
choices of the displacement field D(s). The band pro- 
file can be regarded as either a function of k = k n —k mi 
or equivalently of lu = (E n —E m )/h, related via tu = vk. 
The three band profiles differ in their peak structure, but 
also in their uj — > limits: notably for dilations |M nm | 2 
vanishes in this limit. Our aim is to understand the over- 
all uj dependence, and the small u> behavior in particular. 
For the calculation of band profile we used all 451 eigen- 
states of the 2D quarter-stadium (see Fig. 1) lying be- 
tween 398 < k < 402, found using the method of Vergini 
and Saraceno 0]. For this particular chaotic shape a re- 
markably good basis set (size of order L/A B ) of real and 
evanescent plane waves has been devised which al- 
lows the tension error (defined as the boundary integral 
of ip 2 ) to be typically 3 x 10 -11 in our calculation, (max- 
imum 2 x 10 -10 for any state). The resulting errors in ip 
manifest themselves only when |M„ m | 2 reaches its lowest 
reliable value ~ 10~ 10 , visible as bottoming-out in the 
leftmost point of the inset of Fig. 2. 

In order to understand the quantum-mechanical band 
profile, we can first assume that the eigenstates look like 
uncorrelated random waves. A lengthy but straightfor- 
ward calculation |l2j leads to the result 

I*.,. » «!> ^1 ^ ( „. D) ^, (2) 

where the geometric factors for d = 2, 3, • • • are Sl^ = 
27r,47r--- and <| cos(6>)| 3 > = 4/(3tt), 1/4, • • •. If the dis- 
placement field is normalized such that |D| ~ L, then 
we get |M nm | 2 ~ (Ab/L) 6-1 . Note that the above result 
implies that \M nm \ 2 is independent of w. 
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FIG. 2. The band profile in the 2D quarter-stadium at k ss 400 
for three choices of deformation field: dilation (D), a generic de- 
formation (G), and a generic deformation restricted to parallel 
displacement of the stadium upper edge (Gp). G and Gp are 
chosen to be volume-preserving. In each case, the solid line is the 
average \M nrn \ 2 (estimation error 10%) versus u) = v(k n — k m ), 
with v i— > 1, and the dashed line is the semiclassical estimate 
(Eq.^) (estimation error 3%). We normalized G and Gp so that 
they share the same random wave estimate (Eq.^) as D; this is 
shown as a horizontal dotted line. The system-specific peak due to 
'bouncing- ball' orbits is labelled (bb). The inset is a log- log plot 
with average \M n m\ 2 shown as points. 

To go beyond the random- wave estimate (||), we adopt 
a more physically appealing point of view. We include 
a parametric deformation of the billiard shape via the 
Hamiltonian H(r, p; x) = p 2 /2m + V(r — iD(r)), where 
x controls the deformation. Note that the displacement 
field D is regarded as a function of r. The normal dis- 
placement of a wall element is (n-D)x. The position 
of a particle in the vicinity of a wall element is con- 
veniently described by Q = (s,z), where s is a sur- 
face coordinate and z is a perpendicular 'radial' coor- 
dinate. We set V(r) = Vq outside the undeformed 
billiard; later we take the limit Vq — ► oo. We have 
dTL/dx = — [n(s)-D(s)] Vq8(z). The logarithmic deriva- 
tive with respect to z of an eigenfunction on the bound- 
ary is ip(s)/ip(s). For z > the wavefunction ^(r) is a 
decaying exponential. Hence the logarithmic derivative 
of the wavefunction on the boundary should be equal to 
—y/2mVo/fi. Consequently one obtains (dH/dx) nm — 
— [(hk) 2 /m]M nm , Thus the band profile of M nm is equal 
(up to a factor) to the band profile of the perturbation STL 
due to a deformation of the boundary. See also |p|,[7|JT2| . 

We can now use semiclassical considerations p0| . The 
application to the cavity example has been introduced 
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in |12j]. Here we summarize the recipe. First one 
should generate a very long (ergodic) classical trajec- 
tory, and define for it the fluctuating quantity Tit) — 
—dH.(r, p; x) / dx\ x =o, where the time-dependence of T is 
due to the trajectory (r(£),p(t)). Hence 

T{i) = ^ 2 mv cos(0 co i) D eoJ S(t - t col ) (3) 



where i col is the time of a collision, D cal stands for n-D 
at the point of the collision, and wcos(6' col ) is the normal 
component of the particle's collision velocity. If the defor- 
mation is volume-preserving then (Tit)) = 0, otherwise 
it is convenient to subtract the (constant) average value. 
Now one can calculate the correlation function C(r) of 
the fluctuating quantity T(t), and its Fourier transform 
C(oj) = J C(r) ex-p(iojT)dT. The semiclassical estimate 
for the matrix element is 

where A is the mean level spacing. In practice it is con- 
venient, without loss of generality, to work with units 
such that in (^) the time t is measured in units of length, 
and we make the replacements m i— > 1 and v i— > 1. Then 
(0) can be cast into the form (|M„ m | 2 ) « (Ak/2ir)C(K) 
where Afc is the mean level spacing in k. 

Fig. 2 shows the excellent agreement between the ac- 
tual band profile and that predicted by Eq. (||) for generic 
deformations and dilation. Note that there were no fitted 
parameters in this match. In all estimations of C(u>) we 
have used single trajectories of ~ I0 6 consecutive colli- 
sions. 

Understanding the band profile of \M nm \ 2 has now 
been reduced to a matter of finding a classical theory for 
C(oj). If we assume that Eq.(||) is a train of uncorrelated 
impulses, then its power spectrum would be that of white 
noise, namely C(u>) ~ const. A straightforward calcula- 
tion jl2) then leads to the random wave result (||) already 
presented. However, in reality there are correlations in 
this train, and therefore we should expect C(u>) to have 
some structure on the frequency scale u) ~ l/r col . Look- 
ing at Fig. 2 we see that the white noise expectation is 
reasonably satisfied for one of the 'generic' deformations 
(G), but not in the other two cases (D, Gp). We also see 
non-universal peaks at u> ~ 1/t oo1 ~ 1. We now explain 
that for u> <C 1/ r col there is total failure of the white noise 
result for dilations, as well as for translations and rota- 
tions, and discuss further complications that may arise if 
the billiard system is not strongly chaotic. 
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FIG. 3. The classical power spectrum C(u)) for T(t) correspond- 
ing to a generic deformation (G), dilation (D), translation (T), 
and rotation (R), in the case of the generalized Sinai billiard with 
m = v = 1. Estimation error is 13% for G and R, 20% for D and 
T. The two dotted lines show ui 2 and a; 4 frequency dependence, for 
purposes of comparison. 

In Fig. 3 we display C(ui) for a different billiard shape, 
a generalized Sinai billiard (Fig. 1), chosen because it 
does not suffer from the non-generic marginally stable 
orbits found in the quarter-stadium. Here we see very 
convincing evidence that for small frequencies we have 
C{uo) w const for generic deformation, while C(u>) oc w 4 
for dilation and translation and C(u>) oc w 2 for rotation. 
Thus the white noise expectation is indeed satisfied in the 
u> <C 1/t oo1 regime for generic deformations, but fails for 
dilations, translations and rotations, for which C(lj) — > 
as lj — > 0. This property is known (in the context of 
eigenvalue spectra) as 'rigidity' ||. It implies that the 
train of impulses is strongly correlated, a result which 
at first sight seems inconsistent with the assumption of 
chaotic motion. We will explain that there is no incon- 
sistency here. 

The quantity T(t) = —dH/dx is related to p = 
—dH/dr = — VV, the instantaneous force on the par- 
ticle, by T{t) = — D(r) • p. For translations we have 
D = e, where e is a constant vector that defines a 
direction in space. We can write T{t) — (d/dt) 2 Q(t) 
where Q(t) = —me ■ r. A similar relation holds for 
dilation D = r with Q(t) — —\mv 2 . It follows that 
C(uo) — oj 4 Cg(u), where Cg(uj) is the power spectrum of 
G(t). Assuming that G(t), unlike T(t), is a generic fluc- 
tuating quantity that looks like white noise, it follows 
that C(uj) is generically characterized by u> A behavior for 
either translations or dilations. For rotations we have 
D = e x r, and we can write Tit) — id/dt)Qit), where 
Q{t) = — e ■ (r x p), is a projection of the particle's angu- 
lar momentum vector. Consequently C(w) = oj 2 Cg(o;), 
and we expect C(a>) to be generically characterized by 
(J 2 behavior in the case of rotations. 
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In the previous paragraph we have assumed that 
generic fluctuating quantities such as r 2 and e • r and 
e • (r x p), as well as J- it) for any generic deformations, 
have a white noise power spectrum as u — > 0. Obviously, 
this 'white noise assumption' should be verified for any 
particular example. If the motion is not strongly chaotic, 
meaning that C(r) decays like a power law (say 1/t 1-7 
with < 7 < 1) rather than an exponential, then the uni- 
versal behavior is modified: we may have cj~ 7 behavior 
for small frequencies. For a generic system, for instance 
the generalized Sinai billiard, we do not have this com- 
plication. The stadium example on the other hand is 
non-generic: the trajectory can remain in the marginally 
stable 'bouncing ball' orbit (between the top and bot- 
tom edges) for long times, with a probability scaling as a 
power law in time. Depending on the choice of D(r) this 
may manifest itself in C(r). For example, in Fig. 2 the 
deformation Gp involves a parallel displacement of the 
upper edge, and the resulting sensitivity to the bouncing 
ball orbit leads to large enhancement of the fluctuations 
intensity C{u!—Q), and is suggestive of singular w~ 7 be- 
havior for small u>. 

Finally, consider the time-dependent problem which is 
described by the Hamiltonian H(r, p; x(t)). It is well 
known that under quite general circumstances the dis- 
sipation is ohmic (oc x 2 ). See ]Tl] , |T2| | and references 
therein. If x{t) — Asin(cji), linear response theory gives 
the long-time heating rate d(H)/dt = [i ■ |(uM) 2 . The 
dissipation coefficient /x is determined by the matrix ele- 
ments of ([|), [which up to a factor equals |M nm | 2 ], and 
therefore /i is proportional to C(u>). Our results im- 
ply that fi vanishes in the limit u> — > for translations. 
One should not be surprised Q, since this follows from 
Galilean invariance: One can view the limit u> — > as 
corresponding to the special case of constant x. For con- 
stant nonzero x the particle(s) in the cavity accommo- 
date their motion to the reference frame of the cavity, 
and there is no dissipation. A similar argument holds 
for rotations. On the other hand it is somewhat sur- 
prising that the same conclusion holds for dilations (the 
only other shape-preserving deformation) as well. This 
observation, as far as we know, has not been introduced 
previously in the literature. 

Appendix: There exist a couple of lengthy vector- 
identity proofs ||[l3| of the normalization M nn = 1 for 
the dilation case D = r, for d — 2. Here we present a 
physically illuminating alternative that works for arbi- 
trary d. We use a phase-space-preserving definition of 
dilation operator Uia) = exp(iaG/?i). It is generated by 
the hermitian operator G = | (r • p + p • r) . Applying this 
dilation on wavefunctions gives the expansion: 

U(a)i>(r) « ip(r) + a((d/2)ip + r • Vip) + C(a 2 ) (5) 

The operator also has the effect WrU = e Q r and 
U^pU = c~"p. Consider now any Hamiltonian Ho = 



p 2 /(2m) + V(r). Defining the parameter-dependent ver- 



sion W(r, p;a) = 
ward to obtain 

on 

da 



U(a)Ho(r,p)U(ay , it is straightfor- 



a=0 



= i-__ r . W , 

m 



(6) 



whose matrix elements in the case of the billiard poten- 
tial are idU/da) nm = Hhk) 2 /m) [S nm - M nm }. Thus 
the non-diagonal terms are the same as those of the de- 
formation D = r. The diagonal elements can be calcu- 
lated directly by taking the limit a— >0 of the expression 
i(Uip\n \Uip) - {ip\Ho\ip}) /a. Using (§ and the fact 
that (V>| r -V = — d/2 one can easily show that the 
result equals zero. From here it follows that M nn = 1. 
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